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We analyze the spectral energy distributions (SEDs) of Lyman break galaxies 
(LBGs) at z ~ 1-3 selected using the Hubble Space Telescope (HST) Wide Field 
Camera 3 (WFC3) UVIS channel filters. These HST /WFC3 observations cover 
about 50 arcmin 2 in the GOODS-South field as a part of the WFC3 Early Release 
Science program. These LBGs at z ~ 1-3 are selected using dropout selection 
criteria similar to high redshift LBGs. The deep multi-band photometry in this 
field is used to identify best-fit SED models, from which we infer the following 
results: (1) the photometric redshift estimate of these dropout selected LBGs is 
accurate to within few percent; (2) the UV spectral slope (3 is redder than at high 
redshift (z > 3), where LBGs are less dusty; (3) on average, LBGs at z ~ 1-3 
are massive, dustier and more highly star-forming, compared to LBGs at higher 
redshifts with similar luminosities (0.1L* < L < 2.5L*), though their median val- 
ues are similar within la uncertainties. This could imply that identical dropout 
selection technique, at all redshifts, find physically similar galaxies; and (4) the 
stellar masses of these LBGs are directly proportional to their UV luminosities 
with a logarithmic slope of ~0.46, and star-formation rates are proportional to 
their stellar masses with a logarithmic slope of ~0.90. These relations hold true 
- within luminosities probed in this study — for LBGs from z ~ 1.5 to 5. The 
star-forming galaxies selected using other color-based techniques show similar 
correlations at z ~ 2, but to avoid any selection biases, and for direct compar- 
ison with LBGs at z > 3, a true Lyman break selection at z ~ 2 is essential. 
The future HST UV surveys, both wider and deeper, covering a large luminosity 
range are important to better understand LBG properties, and their evolution. 

Subject headings: galaxies: high redshift — galaxies: fundamental parameters - 
ultraviolet: galaxies — galaxies: evolution 



Introduction 



The high redshift frontier has moved to z > 7 as a result of the high resolution 
near-infrared (NIR) images from the Hubble Space Telescope (HST) Wide Field Camera 
3 (WFC3), and the Lyman break 'dropout' technique. The Ly man break technique was 
first applied to select L yman break galaxies (LBGs) at z ~ 3 ( IGuhathakurta et al.l Il990l ; 
Steidel et al.lll996l . Il999h , and since t hen it has been extensively used to select and study LBG 



cand i dates at redshifts z ~ 3-8 (e.g., 



Bouwens et al.l2007t lHathi etal.l 2008b: Reddv fc Steidel 



20091 : iFinkelstein et al.l I201CH : lYan et al.l 120101 ) . This dropout technique has generated large 
samples of faint star-forming galaxy candidates at z ~ 3-8. However, at highest redshifts 
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(z > 3), it is very difficult to understand the details of their stellar populations using current 
space and ground-based telescopes. Their faint magnitudes make it extremely difficult to 
do spectroscopic studies, and limited high resolution rest-frame optical photometry make 
it challenging to investigate their spectral energy distributions (SEDs). These limitations 
make it imperative to identify and study LBGs at lower redshifts (z < 3). The primary 
reason for the lack of dropout selected LBGs at z ~ 1-3 is that we need highly sensitive 
space-based cameras to observe the mid- to near-ultraviolet (UV) wavelengths required to 
identify Lyman break at z ~ 1-3. The peak epoch of global star-formation rate at z ~ 1-3 
i s now accessible using th e dropout te c hniqu e with the WFC3 UVIS channel 



Hat hi et al 



(bold hereafter H10) and lOesch et all (120 lOf ) have used the HST WFC3 with its superior 
sensitivity to photometrically identify lower redshift (z ~ 1-3) LBGs. Understanding the 
LBGs at z < 3 is vital for two main reasons. First, we need to study the star-formation 
properties of these LBGs, because they are at redshifts corresponding to the peak epoch of 



the global star-form ation rate (e.g.. lReddv et al.ll2008l: iReddv fc Steidelll200a iLv etal.112009 



Bouwens et al.ll2010l ). Second, they are likely lower redshift counterparts of the high redshift 
LBGs — because of their identical dropout selection and similar physical properties — whose 



Labbe et al. 


2010; 


Stark et al. 


2010) 



There are primarily t hree techniques to select star-forming galaxies at z ~ 2: (1) sBzK 
using the B z, K bands, baddi et al.lbooi I200A. (2) BX/BM (using the U, G, R bands, 



Steidel et all 120041 : lAdelberger et al 



2004T) . and ( 3) LBG (using the bands which bracket the 



redshifted Lyman limit, H10; lOesch et al.l |2010| ). All these approaches select star-forming 
galaxies, and yield insight into the star-forming properties of these galaxies, but they have 



differing selection bia ses, and so these samples don't completely overlap (see ILv et al.l 12011 



Haberzettl et al.ll2012l . for details). Therefore, it is essential to apply identical selection crite- 
ria at all redshifts to properly compare galaxy samples and accurately trace their evolution. 
The LBG selection is widely used to select high redshift (z > 3) galaxies, and to do equal 
comparison with these galaxies, here we investigate physical properties of LBGs at z < 3. 

H10 used UV observations of the WFC3 Science Oversight Committee Early Release 
Science extragalactic program (PID: 11359, PI: O'Connell; hereafter "ERS"), which covers 
approximately 50 arcmin 2 in the norther n-most part of the Great Observatories Origins Deep 
Survey (GOODS; iGiavalisco etaDl2004h South field, to identify LBGs at z ~ 1-3. The high 
sensitivity of the WFC3 UVIS channel data ( jWindhorst et al.l 1201 ll ). along with existing 
deep optical data obtained with the Advanced Camera for Surveys (ACS) as part of the 
GOODS program are ideal to apply dropout technique in observed UV filters to select LBG 
candidates at z ~ 1-3. In this paper, we use this H10 sample of LBGs to investigate their 
physical properties by fitting stellar synthesis models to their observed SEDs. 
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This paper is organized as follows: In § (2J we summarize the WFC3 ERS observations, 
and discuss our LBG sample at z ~ 1-3 as well as the comparison sample of LBGs at z ~ 4-5. 
In § [21 we fit observed SEDs of LBGs at z ~ 1-3 and z ~ 4-5 to stellar population synthesis 
models, and discuss the best-fit parameters (redshift, UV spectral slope, stellar mass, stellar 
age, and star-formation rates) obtained from these SED fits. In § HI we discuss correlations 
between best-fit physical parameters and their implications on our understanding of LBGs. 
In § we conclude with a summary of our results. 

In the remaining sections of this paper we refer to the HST /WFC3 F225W, F275W, 
F336W, F098M, F125W, F160W, filters as U 225 , U 275 , U 336 , Y 098 , J 125 , H 160 , to the HST /ACS 
F435W, F606W, F775W, F850LP filters as £435, ^eoe, ^775, z 850 , and to the Spitzer /IRAC 
3.6 fim, 4.5 /im filters as [3.6], [4.5], respectively, for convenience. We assume a Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) cosmology with Q m =0 .27A, £^=0-726 and i fo=70.5 km 
s _1 Mpc -1 , in accord with the 5 year WMAP estimates of iKomatsu et al.l (120091 ). This cor- 
responds to a look-b ack time of 10.4 Gyr at z ~ 2. Magnitudes are given in the AB^ system 
Joke fc Gunnlll983l l 



2. Observations and Sample Selection 



The WFC3 ERS observations (IWindhorst et al.ll201ll ) were done in both the UVIS (with 
a FOV of 7.30 arcmin 2 ) and the IR (with a FOV of 4.65 arcmin 2 ) channels. Here, we briefly 
summarize the UV imaging observations. The WFC3 ERS UV observations were carried 
out in three broad-band filters ^225, ^275 and U 336 . The U 225 and U 275 filters were observed 
for 2 orbits (~5688s) per pointing, while the U 33 q filter was observed for 1 orbit (~2778 
s) per pointing, for a total of 40 orbits over the full ERS field (8 pointings). We used the 
existing GOODS v2.C0 reduction of the ACS images in four optical bands (-B435, V606, 2775, 
z 85o), which were re-binned to a pixel size of 0.09". To match the ERS IR (logs, ^125, #160) 
and re-pixellated ACS optical images, the UV mosaics have a pixel scale of 0.090" pix^ 1 
and cover ~50 arcmin 2 area of the GOODS-South field. Details of these observations and 



reduction process are described in I Windhorst et al.l (120111 ). 



The combination of the three WFC3 UV filters a nd the four ACS o ptical filters provide 



an excellent ability to select LBGs at z ~ 1-3 (H10; lOesch et al 



2010), using the dropout 



technique to detect the Lyman-break at rest-frame 912 A (lMadaulll995l ). H10 used dropout 
color selection technique in the ERS UV field to identify three sets of UV-dropouts — t/225- 
dropouts, t/275-dropouts and [/336-dropouts — which are LBG candidates at z ~ 1.6, 2.2 and 



1 http://archive.stsci.edu/pub/hlsp/goods/v2/ 
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2.6, respectively. They found 66 U225-, 151 ^275- and 256 l^-dropouts to a magnitude limit 
of AB ~ 26.5 mag, respectively. 

In this paper, we start with the H10 LBG sample. Our goal is to investigate SEDs 
of reliable LBG candidates with at least 10-band HST coverage (augmented by additional 
data as described in § [3]) from the H10 sample, so we apply the following filtering criteria. 
First criterion is the availability of the WFC3 IR (I098) J125, -H160) data. H10 used the 
WFC3 UVIS and ACS data to select LBGs at z ~ 1-3. The WFC3 UVIS channel has a 
larger final ERS mosaic than the WFC3 IR channel, so we exclude LBG candidates that 
don't have WFC3 IR data from our SED analysis. This criterion reduces the H10 sample 
size by about 10%. Secondly, galaxies with poor SED fits (see § as measured by their 
larger x 2 were excluded from the sample. The galaxies which fail the SED fit are usually 
fainter, does not have all NIR photometric data or have highly uncertain NIR photometry, 
has poor best-fit SED (indicated by high x 2 ) and could have a primary lower-redshift (z < 1) 
solution. This criterion removes additional ~10% of galaxies from the H10 sample. This 
fraction of catastrophic x 2 outliers is consistent with the outlier fraction in the photometric 
redshift distribution of the dropouts in the H10 sample. The final sample of LBGs for the 
SED analysis — after applying above mentioned criteria — is 47 U225-, 126 U275- and 213 
C/ 3 3 6 -dropouts. 

To compare SED properties of LBGs at z ~ 1-3, we select -B435- and V<306-dropouts in 
the WFC3 ERS field. These dr opouts — LBG candi dates at z ~ 3.7 and 4.7, respectively 



were selected following the iBouwens et al.l (120071 ) selection criteria. The goal of this 



paper is to compare the HST/WFC3 selected UV-dropout galaxies at z — 1-3 with similar 
galaxies at higher redshifts (z ~ 4-5). The limited area and depth of the ERS data puts 
brightness limitations on our UV-dropout selection, which was confined to comparatively 
brighter part of the rest-UV luminosity function (around the knee and brighter as shown in 
H10). Therefore, we have selected i?4 35 -dropouts and V^-dropouts in the ERS field whose 
luminosities are similar to the UV-dropouts. The luminosity range is 0.1L* < L < 2.5L* 
(based on L* corresponding to M=-21 mag). Applying similar filtering criteria as LBGs at 
z ~ 1-3, we have a comparison samp le of 155 -B435- drop outs and 27 VeQ6-dropouts in the 



ERS field. Based on IXue et al.l (120111 ) X-ray catalog, there are four active galactic nuclei 
(AGN) in the H10 LBG sample (z < 3), three AGN in the 5 435 -dropout sample, and none 
in the V(306-dropout sample. These small number of X-ray AGN does not affect our results 
or conclusions. 

All subsequent analysis in this paper is done identically on these 5 samples (U225-, U275-, 
U336-, -B435-, and V(306-dropouts) for proper comparison. To show general evolutionary trends, 
we combine three dropout samples from H10 as a UV-dropout sample (z ~ 1-3), and two 
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high redshift samples as a I?435-,V606-dropout sample (z ~ 4-5). 



3. Spectral Energy Distributions 



The Le PHARE software package (jArnouts et al.lll999c lllbert et al.ll2006l ) was used to mea- 
sure the photometric redshifts, and to fit the broadband SEDs of LBGs. The primary goal 
of SED fitting is to find the best-fitting synthetic stellar population model to the observed 
photometry. From this best-fit model, we can estimate the redshift, stellar age, stellar mass, 
star-formation rate (SFR), dust extinction, and ot h er ph ysical properties of each galaxy. 
We use the 2007 version of the iBruzual &: Charlotl (120031 . hereafter CB07) models, which 
has improved prescription of thermally pulsating AGB stars. We generated a set of stel- 
lar population models assuming a Salpeter initial mass function, and varying the redshift 
(z = 0.1-6.0, <5z=0.1, though a parabolic interpolation is used to refine the photometric 
redshift solution within Sz intervals), metallicity (0.2, 0.4 a nd 1 Zm), age (1 My r < t < tn), 
dust extinction (0 < E(B-V) < 0.7 mag, using a modified ICalzetti et al.l l2000l attenuation 
law), and e-folding timescale (y=0. 1,0. 3, 1,2, 3, 5, 10, 1 5, 30 Gyr) for a star-formation history 
(SFH) ocexp(-t/r). The Le PHARE code assumes the iMadaul (119951 ) prescription to estimate 
inter-galactic medium (IGM) opacity. The model that gives the lowest x 2 is chosen as the 
best-fit SED. 

The contribution of major emission lines in different filters can be included in the models 
using the Le PHARE code. Neglecting emission lines during the SED fit ting; process can over- 
estimate the best-fit stellar ages and masses b y as much as 0.3 dex (e.g. JSchaerer fc de Barros 
20091 ; iFinkelstein et al.ll201l[|Atek et al.lboill ). The Le P HARE code accoun ts for the contribu- 
tion of emission lines with a simple recipe based on the iKennicuttl ( 119981 ) relations between 
the SFR and UV luminosity, Ha and [Oil] lines. The code includes the Lya, Ha, H/3, 
[Oil], QUI [4 959] and OI I I[500 7] lines with different line ratios with respect to [Oil] line, as 
described in lllbert et all (I2009I ). 



The observed photometry is available in up to 13 filters: three HST /WFC3 UVIS, four 
HST/ ACS, three JZST/WFC3 IR, one VLT Ks, and two Spitzer /IRAC [3.6], [4.5] bands. 
We perform matched aperture photometry in 10 HST bands as discussed in H10, while 
we use VLT and S pitzer photometry from the publicly available GOODS-MUSIC catalog 
( jSantini et al.ll2009l ). The photometry in MUSIC catalog has accurate PSF-matching of space 
and ground-based images of different resolution and depth. Figure [1] shows example best-fit 
SEDs for LBGs at z ~ 1-3, and the comparison sample of LBGs at z ~ 4-5. 
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3.1. Photometric Redshifts 



One of the free parameters during the SED fitting process is the redshift. To assess the 
accuracy of our SED-based photometric redshifts (z p h) at z ~ 1-3, we compare them with 
the spectr oscopic redshifts (z. 1P ) from various VET/Magellan campaigns in the GOODS-S 



field fe.g.. iGrazian et al 



its \Zxp) irom various y Ll / Magellan campaigns m tne LrU> 
200d; iRavikumar et al.l 120071; IVanzella et al]|2008l : IWuvts et al 



2008 



Balestra etallboiol : ICooper etallboij ). We find that only a small number (< 30%) of H10 



dropout sample has spectroscopic redshifts, most likely due to the lack of strong features 
in observed 4500-9000 A range at 1 < z < 3, where most ground-based spectrographs on 
large telescopes are optimized. We matched 91 spectroscopic redshifts for the UV-dropout 
sample [z ~ 1-3) selected with the criteria discussed in §[2J Figure [2] shows the comparison 
between the SED based photometric redshifts and the publicly available spectroscopic red- 
shifts. The catastrophic outliers — shown by concentric circles in Figure [2] — have quality 
flags that indicate the spectroscopic redshift is unreliable (in most catalogs C or worst). So 
it is likely that these spectroscopic redshifts are not correct and hence, redshift comparison 
for these objects is not credible. The histogram in Figure [2] shows the distribution of photo- 
metric redshift uncertainties 5z=(z sp — z p h/l+z sp ). Based on this distribution, we estimate 
cr(fc)~0.05, and < 5z >=-0.03. The fraction of catastrophic outliers (>3<r) is ~7%, exclud- 
ing objects with u nreliable spectroscop ic redshifts. Our photometric redshift uncertainties 
are consistent with iDahlen et al.l (120101 ) , who used the deepest a nd the most comprehensive 
photometric data in the GOODS-S field. iHaberzettl et al.l (120121 ) selected fairly bright LBGs 
at z ~ 2 using the GALEX data, and found similar photometric redshift uncertainties and 
outlier fraction for their dropout sample. The photometric redshift uncertainty in the im- 
plied redshift is also consistent with the dropout selection method applied to select these 
galaxies. The dropout selection technique uses the location of a spectral break within a 
photometric bandpass (filter), and therefore, the redshift uncertainty depends on the width 
of the bandpass, and could be as high as ~0.5 in z. 

The distribution of photometric redshift uncertainties (5z) is asymmetric, even after 
excluding objects with unreliable spectroscopic redshifts. There are more galaxies in the 
distribution with spectroscopic redshifts lower than their photometric redshifts i.e., (z sp — 
z p h < 0). Detail investigation of each ground-based spectra (if available) is needed to figure 
out what is causing spectroscopic redshift to be lower than photometric redshift. Such an 
investigation is beyond the scope of this paper, but we should p oint out that such asym metric 



distribution is also observed for GALEX-selected LBGs (e.g.. IHaberzettl et al.ll2012l ). and is 



totally consistent within the estimated photometric redshift uncertainties. 
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3.2. UV Spectral Slope 8 



The UV spectral slope 8 is determined from a power-law fit to the UV continuum 
spectrum ( ICalzetti et al.l Il994h . fx oc X 13 , where f\ is the flux density per unit wavelength 
(ergs s _1 cm" " 2 A" 1 ). We use the best-fit SEDs of dropout selected LBGs to estimate their UV 
spectral slope j3 by fitting a straight line between rest-frame 1300 and 1900 A in their model 
spectrum. Th i s wav elength range covers 7 out of 10 spectral fitting windows identified by 
Calzetti et al.l ( 119941 ) to estimate the UV spectral slope. This wavelength range is also ideal 
for comparing 8 values at higher redshifts, because those are usually measured between rest- 
frame 1600 and 2000 A. Figure [3] shows the slope- fitting method applied to the best-fit SEDs 
to estimate 8, where the solid line is the best-fit UV spectral slope, the dashed line is the 
best-fit SED, and the black filled circles are observed magnitudes. By selection (§|2J), we only 
consider galaxies with good SED fits so the choice of model should not affect the (3 estimate. 
Though uncertainties in observed photometry could affect the best-fit SED parameters, and 
hence, the (3 estimate. In Figure [3] , we also quote typical intrinsic uncertainty in /3 for 
galaxies at different redshifts. We estimate (3 for each galaxy, and then fit a Gaussian to the 
8 distribution to find median (and sigma) value in each redshift bin. Tabled] shows median 8 
values and their corresponding uncertainties for the UV-dropout, and the _B 43 5-,V606-dropout 
samples. 

The evolution in the UV spectral slope 8 as a function of redshift may indicate change 
in stellar populations of galaxies over cosmic ti me. We compare our 8 values with the highe r 



redshift measurements from the literature (e.g.. iBouwens et al.ll2012l ; iFinkelstein et al.ll2012l ). 
Figure 0] shows the UV spectral slope 8 as a function of redshift. Blue filled squares are 
median 8 values measured between rest-frame 1300 and 1900 A for our dropout samples. 
To test how 8 measurements are affected by the selection of rest-frame UV wavelength 
range, we also measured (3 between rest-frame 1300 and 3400 A, which are shown by blue 
open squares in Figure HI Both 8 values agree within la u ncertainties. The red filled 
diam onds are measurements from IBouwens et al.l ((2009, 201 2J) an d purp l e fille d circles are 
from Finkelstein et al. ( 2012 ). The 8 values from Bouwens et al. ( 20091, 2~012h are f or the 
galaxies around M* v , which is consistent with our sample, while IFinkelstein et al.l (120121 ) 
measurements are based on all galaxies extending to those fainter than M* v in their respective 
redshift bins. The uncertainties on the median values of (3 are the stan dard error of th e mean 



in th e case of IFinkelstein et al.l (120121 ) and our measurements, while IBouwens et al.l (12009 



20121 ) uncertainties represent la scatter. For comparison, our estimated la scatter in median 
8 values are listed in Table [TJ 

Figure H] shows that the median values of 8 decreases as redshift increases (/3 ~-1.6 
at z ~ 1.6 to j3 ~— 2.4 at z ~ 8), which could imply variations in one or more physical 
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properties of LBGs as a function of redshift. Figure [5] shows evolution in (3 as a function 
of best-fit SED parameters (stellar mass, stellar age, dust content, SFR) and redshift. The 
lowest redshift bin (z — 1.6) is shown by the smallest circles, and the highest redshift bin 
(z — 4-5) is shown by the largest circles. The largest change (factor of ~2 or 0.3 dex) is 
seen in the dust content E(B-V) of galaxies as f3 changes from -1.6 (at z ~ 1.6) to —1.9 (at 
z ~ 4-5), while other three parameters vary much less than a factor of 2. This could imply 
that change in the dust content of LBGs has largest effect on the UV spectral slope (3, and 
any variation in (3 as a function of redshift could most likely be due to changing dust content 
of galaxies. Therefore, based on our estimates, as shown in Figure HJ we could say that 
LBGs at lower redshift (z — 1.6) have more dust than LBGs at higher redshift {z — 5). This 
trend of j3 is consistent with pre vious studies, which argued that galaxies at z ~ 6 te nd to be 



bluer than those at z ~ 3 (e.g.. IStanway et al.ll2005t iBouwens etalj|2006j ; lHathi et aLj|2008a 



Wilkins et al.l 1201 ll ). Those (3 measurements, on uniformly selected LBGs, were limited to 
LBGs at z > 3, and our results below z ~ 3 extends this observed trend to z ~ 1.5. 

The evolution in the UV spectral slope could also be due to changing star-formation 
history, initial mass function (IMF), and/or metallicity. These effects are believed to be 
much smaller than the effects from changing dust content of the galaxy. Many authors 



have i nvestigated various stellar population models to estimate these effects. IBouwens et al. 



(120121 ) explored sensitivity of the UV-continuum slope f3 to changes in the mean metallicity, 
age, or dust extinction by choosing one fiducial model as a benchmark, and then changing 
various model parameters to assess changes in (3. They conclude that a factor of 2 (or 0.3 dex) 
changes in metallicity, age or dust content result in 0.07, 0.15, 0.35 changes in the UV spectral 
slope (3, respectively. This implies that changes in the dust content have much larger effect 
on the UV-continuum slop e than similarly-sized changes in the age , metallicity, or the s tellar 
IMF. Similar studies (e.g., iLeitherer et al.l [l999l : lHathi et aDl2008al : Iwilkins et allboilh have 
come to the same conclusion that though (3 could be affected by various stellar population 
properties, the change in dust content of galaxies is the predominant effect which causes f3 
to change. We should also emphasize that it is very challenging to completely understand 
these various effects based on observations only, rigorous modeling and/or simulations are 
required to fully assess the contributions of these various effects on the UV spectral slope f3. 



3.3. Stellar Population Properties 

We compare observed SED with a suite of model templates from CB07 to find the best- 
fit model through \ 2 minimization. All SED parameters are fit simultaneously. The best-fit 
model allows us to estimate photometric redshift (as shown in § 13.11) . and physical properties 
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of stellar populations such as stellar age, stellar mass, dust extinction E(B-V) and SFRs for 
each galaxy. The estimated uncertainties (~0.3-0.4 dex) in stellar ages, masses and SFRs 
are estimated by marginalizing the uncertainties in observed photometry and redshift. 

One of the main limitations of SED fitting is the need to assume a SFH, which cannot 
be reliably constrained for a galaxy from limited photometric data points. We have assumed 
an exponentially declining SFH. Different SFHs (e.g., rising, constant, declining) introduce 
systematic u ncertainties in the stellar mass determ inations, mostly at redshift greater than 



leg 



dex 



e.g., 



Lee et al. 



20ld 



Finlator et al. 



Papovich et al.l 1201 ll ). These uncertainties are typically <0.3 



20071 ). and are within our estimated uncertainties. Stellar ages are 
highly sensitive to the assumed SFH. Any prior star-forming episode can be overshadowed by 
newly born stars from the most recent star-formation, totally neglecting possible existence 
of older population in a given galaxy. Therefore, based on assumed histories, it is possible 
to get an older or a younger age for the same galaxy. Hence, interpreting the stellar ages 
derived from SED fitting can be tricky, and at the very least, the uncertainties on the stellar 
ages could be much larger than estimated uncertainties (~0.3-0.4 dex). Because of these 
issues, in subsequent analysis, we will not elaborate on stellar population ages and focus on 
other physical properties. 

Figure [6] shows the distributions of stellar age, stellar mass, SFRs, and E(B-V) for LBGs 
at z ~ 1-3 (black), and the comparison sample of LBGs at z ~ 4-5 (red). The median values 
are shown by dashed vertical lines and la uncertainties in these distributions for LBGs at 
z ~ 1-3 are shown by an error bar at the top of the black histogram. A two-sided K-S test - 
in each panel — indicates a probability less than 0.006 that the distributions (red and black 
histograms) are drawn from the same parent distribution. Figure [6] shows a general trend 
that — on average — higher redshift LBGs have low SFRs, less dust, and are less massive 
than their lower redshift counterparts, though median values of two distributions (red and 
black) are similar within la uncertainties. This r esult is in good ag reement with previous 



studies comparing LBGs at z ~ 3 and z ~ 5 (e.g., IVerma et al.l 120071 ). The average E(B-V) 
at z ~ 2 is consiste nt with studies based on star -forming; galaxies selected using BX/BM 
color technique (e.g.. lErb et al.ll2006l ; ISawickil 120121 ). The distribution of E(B-V) completely 
agrees with the UV spectral slope evolution as discussed in § 13.21 implying that the LBGs 
at z ~ 1-3 are more dusty (redder) compared to LBGs at z ~ 4-5. 
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Results and Discussion 



4.1. Stellar Mass vs UV Luminosity Relation 



The rest-frame UV light traces recent or instantaneous SFR, while rest-frame optical 
and NIR data help us to estimate stellar masses of galaxies. If the galaxy stellar mass and 
UV luminosity are related then we can directly use rest-frame UV light to estimate stellar 
mass without needing rest-frame optical/NIR data. Figure [7] shows stellar mass of LBGs at 
z ~ 1.5-5 as a function of their UV absolute magnitude. These quantities are based on best- 
fit SEDs, and their typical uncertainties are shown in the lower-left corner. The dotted lines 
are best-fit l ine ob tained by keeping the logarithmic slope fixed at 0.46, which was estimated 
by ISawickil ( 120121 ) for star- forming galaxies at z ~ 2. The dot-dash lines show the scatter 
from the best-fit line, which is ~0.3 dex for LBGs at z ~ 1-3 and about 0.2 dex for LBGs 
at z ~ 4-5. We also tested the validity of this relation by fitting the slope of the line rather 
than fixing it. We find that the fitted slope is in the range of 0.42±0.06 for our LBG samples, 
which is consistent with 0.46 within the estimated la scatter in this relation. Therefore, we 
find that a proportionality relation between these two parameters with a logarithmic slope 
of 0.46 provides a good fit to the data. The stellar mass es of the brighter L BGs — with 



UV luminosities near the L* value of LBGs at z ~ 3 from Steideletal 



(Il999h — are about 



a factor of 2 lower than 10 10 M estimated by iPapovich et al.l (j200ll ). This discrepancy, 
though within our estimated uncertainties, could be due to the fact that we include emission 
lines in our SED fitting which could affect stellar masses by as much as a factor of ~2. The 
stellar mass-U V luminosity relation is fairly tight with a small scatter (< 0.3 dex) , which is 
consistent with other studies at similar redshifts (e.g., IPapovich et al.l 1200 ll ). and it points 
to a nearly constant mass-to-light ratio (log(M/L) ~-0.5) for LBGs between z ~ 1.5 and 
5. The tightness/lower scatter of the stellar mass-UV luminosity relation in Figure [7] could 
be — in part — due to the fact that both these quantities are output parameters from the 
best-fit SEDs, and therefore, it is possible that these parameters are not totally independent. 
We have addressed this issue and discussed its implications in § 14.31 A similar correlation 
between stellar m ass and absolute magnitude has been reported for LBGs at z ~ 5-6 by 
Stark et all (|2009h . 



Figure [7J shows that LBGs at z ~ 1.5-5 follow similar linear correlation between stellar 
mass and UV absolute magn itude (within uncerta inties) for M uv between -19 and -22.5 mag. 
It is important to note that IShapley et al.l (120051 ) does not find any correlation between the 
stellar mass and UV absolute magnitude for star-forming galaxies at z ~ 2 with stellar masses 
> 10 10 VLry This co uld be due to different color-selection technique (BX/BM) used by the 



Shapley et al.l (120051 ) to select star-forming galaxies, whose physica l properties could diffe r 
from the dropout selected LBGs at these masses (e.g.. lLy et al.ll2011uHaberzettl et al.ll2012l ). 
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It is also possible that their sample — which consists of spectroscopically confirmed bright 
galaxies with stellar masses greater than or equal to 10 10 M — has more massive galaxies 
than our sample and it is uncertain how massive galaxies would follow this correlation. The 
ERS observations are too limited in area and depth to cover a larger luminosity range, so 
we cannot predict how this relation will evolve for luminous (M uv <-22.5 mag) or dwarf 
(M uv > -19 mag) galaxies at these redshifts. 



4.2. SFR vs Stellar Mass 



The correlation between the current SFR and stellar mass in star-forming galaxies, 
also known as 'main sequence of star- formation ' (MS ), has been observed at z < 2 (e.g., 
Noeske et al.l 120071 ; lElbaz et al.l 120071 ; iDaddi et ajj 120071 ). These studies have shown that the 
MS relation seems to be not evolving strongly with redshift, but the zeropoint does: that is 
high redshift (z ~ 2) star-forming galaxies are forming stars at a higher rate than similar 
mass local galaxies. In Figure El we investigate this relation and star-formation histories 
for LBGs at z ~ 1.5-5. These quantities are based on best-fit SEDs, and their typical 
uncertainties are shown in the lower-right corner. The dotted lines are best-fit line obtained 
by ke eping the logarithmic slope fixed a t 0.90, estimat ed for star-forming galaxies at z < 2 



(e.g.. 



Elbaz et al.ll2007t IDaddi et al.ll2007l ; ISawickill2012l ). The dot-dash lines show the scatter 
from the best-fit line, which is ~0.6 dex for LBGs at z ~ 1-3 and about 0.4 dex for LBGs 
at z ~ 4-5. We also obtained the best-fit logarithmic slope for this relation, and found 
the fitted slope in the range of 0.81±0.30 for our LBG samples, which is consistent with 
0.90 within the estimated la scatter in this relation. We find that a proportionality with 
a logarithmic slope of 0.90 provides a good fit to the data with few outliers at stellar mass 
greater than 10 10 M . 



Finlator et al.l (120061 ) have shown that tight relation exists between SFR and stellar 
mass for galaxies at z ~ 4 using the cosmoloKical hyd r odynamic simulation s, which is also 



consistent with the observations (IBouwens et al.l 120121 ) . iFinlator et al.l ( 120061 ) also point out 
that the scatter in the Figure |8] could be a measure of SFR 'burstiness' as a function of stellar 
mass. This means that the linear relation (with a logarithmic slope of ~0.90) indicate an 
average SFR for a given stellar mass, but galaxies can also experience bursts of up to two 
times the average SFR value at the same stellar mass as shown by the scatter. The scatter 
in the SFR versus stellar m ass relation for L BGs at z ~ 1-3 is slightly larger than ~0.3 
dex — observed at z ~ 2 by IDaddi et al.l (120071 ) — possibly because of few galaxies forming 
a sharp edge towards high SFR values, as seen in the relation for U225- and [/275-dropouts 
(upper panel in Figured]). These galaxies have low stellar ages (less than 10 Myr), which 
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could be highly uncertain as discussed in § 13.31 It is also possible t hat this edge could be an 



artifact due t o lower limits on the model parameters r and t (e.g.. lHainline et al.l 120121 ) . We 
also note that iMcLure et al.l ( 1201 ll ) argue that the tightness in the SFR-stellar mass relation 
depends on the assumed SFH. The scatter in this relation is much less for a constant SFH, 
while it is much larger for other SFHs. Therefore, it is also likely that the larger scatter we 
see in Figure [8] could be due to different SFHs. 

Figure E] shows that, though our data has little more scatter compared to the MS 
relation at z < 2, the majority of our galaxies fall on to this relation char acterized by a 
l ogari thmic slope of 0.90. A similar correlation is observed at z ~ 6-8 by IMcLure et al 



( 1201 ll ). and supported by cosmological hydrodynamic simulations of iFinlator et al.l ( 1201 ll ). 
Our observations confirm this MS relation for star-forming galaxies from z ~ 1.5 to 5, 
implying that — on average — their star-formation histories are similar. 



4.3. Implications 



In previous sections, we have shown that LBGs at z ~ 1-3 — on average — are massive, 
dustier, and have higher star-formation rates than LBGs at z ~ 4-5 with similar luminosi- 
ties, though it should also be not ed that they are not very different within estimated la 
uncertainties. As pointed out by iPapovich et al.l (120111 ). the number densities of galaxies 
at fixed luminosity could change substantially over this redshift range, which could lead to 
potential biases when comparing galaxies at different redshifts. However, the general trends 
we observe in stellar masses, SFRs, and dust extinction are supported by other independent 
means. The characteristics UV luminosity (L* v ) is increasing as a functio n of redshift from 
z ~ 8 to 2 (e.g., H10), which implies increase in SFRs with time, while iFinkelstein et al. 
(120101 ) have shown that stellar masses for M* v LBGs grow from z ~ 8 to 2. The UV spectral 



slope j3 shows evolution as a function of redshift (Figure HJ), which could indicate lower dust 
content at higher redshifts. The higher dus t content in LBGs at lower redshift is also i n 
accordance with the studies at z 



1 (e.g. 



Burgarella et al.l 120071 iBasu-Zvch etal.l |2011[ ). 
while the IVerma et al.l ( 120071 ) supports the lower dust content in LBGs at z ~ 5. Therefore, 
the ensemble properties of LBGs in our sample are in general agreement with the expected 
results. 

The stellar mass-UV luminosity relation (Figure [7]) and the SFR-stellar mass relation 
(Figure |8]) are based on measurements from best-fit SEDs, therefore, it is possible that these 
quantities are not totally independent, which might affect their observed correlations. To 
investigate this, we show distributions of mass-to-light (M/L; Mass/L uv ) ratios and specific 
SFRs (SSFR; SFR/Mass) in Figure M The black (red) histograms show distribution for 
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LBGs at z ~ 1-3 (z ~ 4-5), and the median values are shown by dashed vertical lines. The 
median values of M/L ratio and SSFR for LBGs at z ~ 4-5 are slightly lower than that at 
z ~ 1-3, but are still consistent within the lex uncertainties as shown by the error bar on the 
top of the black histogram. A two-sided K-S test — in each panel — indicates a probability 
less than 0.05 that the distributions (red and black histograms) are drawn from the same 
parent distribution. The constancy of the M/L ratio and SSFR between z ~ 1.5 and 5 agrees 
very well with the constant slope we find in Figure [7] and Figure |8] for our sample of LBGs, 
though with a slightly larger scatter. 



Stellar masses of LBGs at z 



Elbaz et al. 


2007; 


Daddi et al. 


2007; 


Sawicki 


2012) 



1-3 are generally well correlated with UV absolute 
ning galaxies at similar redshifts (e.g., 



mass assemb ly and SFH for t hese galaxies, but the exact nature of SFHs is still not clearly 



understood. iPapovich et all ( 120111 ) showed that the cosmologically averaged SFRs of star- 
forming galaxies at 3 < z < 8 — at constant co-moving number density — increase smoothly 
from z = 8 to 3, and the stellar mass gro wth in these gal axies is consistent with this derived 
SFH. The scenario of rising SFH (see also lLee et aJboioh is also supporte d by recent results 



from the cosmological hydrodynamic simulations (e.g.. lFinlator et al.ll201ll ). The models with 
rising SFHs conflicts with the assumptions that t he SFR in distant galaxies is either con- 



stant or decreasin g exponentially with time (e.g., IPapovich et all 120051 ; IShapley et all 12005 
Labbe et all 120101). Though, we remind the reader that the models with rising SFHs advo- 
cated by IPapovich et al.l (120111 ) and others correspond to a cosmologically averaged SFHs 
for typical galaxies, and not individu al galaxies, because th ey c ould involve ran dom events 
that changes their instantaneous SFR. IPapovich et al.l (120 111 ) and iLee et al.l (120101 ) also argue 
that rising SFHs are most beneficial to higher redshift [z > 3) galaxies. We find that for our 
assumed SED model parameters, the LBGs between redshift z ~ 1.5 and 5 — on average - 
have similar SFHs, though the precise nature of SFHs at all redshift is still under debate, 
and could also affect the SFR-stellar mass correlation. 

Our analysis demonstrates that the dropout selected galaxies at z ~ 1-3 — within lumi- 
nosities probed here — show similar correlations between physical parameters (SFR, stellar 
mass, UV luminosity) as other star-forming galaxies selected u sing diffe r ent c olor criteria 
(e.g., sBzK, BX/BM) at z ~ 2. This is consistent with the iLy et al.l ( 1201 ll ) conclusion 
that majority (~80-90%) of the dropout selected galaxies overlap with other color selected 
star-forming galaxies with stellar masses less than 10 10 M . The stellar mass range for our 
current sample is between ~10 8 and ~10 10 M Q , with sample completeness around 10 9-9 - 5 
M . Significant differences between the dropout selected sample and other color selected 



samp les of star-forming galaxies at z ~ 2 exists for massive galaxies (>10 10 M ; iLv et al 



201ll ). Therefore, it is vital to use uniform selection technique at all redshifts to avoid any 
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selection biases. The Lyman break dropout technique is the most convenient and widely 
used method to select galaxies at z > 3, and we have shown that LBGs at z — 1-3 selected 
using this dropout technique have similar physical properties (within uncertainties) as LBGs 
at z ~ 4-5 with similar luminosities. Hence, LBG selection at z < 3 is important to under- 
stand properties of LBGs and properly investigate their evolution as a function of redshift. 
The validity of LBG properties over wide luminosity and mass range can be in vestigated in 



detai l with the upcoming an d future WFC3 UV surve ys such as CANDEL S (IGrogin et al. 



20 111 ; iKoekemoer et al.ll2011[ ) and the WFC3 UV UDF flRafelski et all 12012 . PI: H. Teplitz). 



5. Summary 

In this paper, we investigated stellar populations of LBGs at z ~ 1-3 selected using 
HST /WFC3 UVIS filters in the GOODS-S field. We used deep multi-wavelength observa- 
tions from the HST, VLT, and Spitzer to compare observed SEDs with the spectral synthesis 
models to infer physical properties (stellar masses, stellar ages, SFRs, and dust extinction) of 
these LBGs. We also compared these LBGs with their higher redshift (z ~ 4-5) counterparts 
with similar luminosities (0.1L* < L < 2.5L*). Our results can be summarized as follows: 

• We obtain reliable (cr(z sp —z ph /l+z sp )~0.05) photometric redshifts for dropout selected 
LBGs at z ~ 1-3 based on 10-13 band SEDs. 

• The UV continuum slope f3 for LBGs at z ~ 1-3 is redder (j3 ~ -1.6 at z ~ 1.6) 
compared to their higher redshift counterparts ((3 ~-2.4 at z ~ 8), implying higher dust 
content in these LBGs. 

• On average, LBGs at z ~ 1-3 are massive, dustier and more highly star-forming com- 
pared to LBGs at z ~ 4-5, though their median values are very similar within estimated la 
uncertainties. This similarity emphasizes the importance of identical Lyman break selection 
technique at all redshifts, which selects physically similar galaxies. 

• The stellar mass-absolute UV magnitude relation for LBGs between z ~ 1.5 and 
5 show linear correlation with a logarithmic slope of ~0.46, while the SFR-stellar mass 
relation show similar correlation with a logarithmic slope of ~0.90. To properly compare 
and interpret such relations at higher (z > 3) redshift, and to avoid any selection biases due 
to different selection techniques, a true Lyman break selection is required at z ~ 2. 

• We need larger HST UV surveys to cover full range in luminosity/mass and better 
understand LBG properties, and their evolution. Both deeper and wider UV surveys are 
needed. The wider one to probe the high mass end, while the deeper one will probe the sub- 
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L* population. A large number of HST orbits have been used for dropout selected galaxies 
at z > 3, and the lower redshift regime needs to be explored in a comparable manner. 

We thank the referee for helpful comments and suggestions that significantly improved 
this paper. This paper is based on Early Release Science observations made by the WFC3 
Scientific Oversight Committee. We are grateful to the Director of the Space Telescope 
Science Institute for awarding Director's Discretionary time for this program. Support for 
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Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. This research was (par- 
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Fig. 1. — [ Top panel] Example best-fit SEDs of LBGs at z ~ 1-3. We have used the Le 
PHARE code (jArnouts et al.lll999t lllbert et al.ll2006l ) to compute photometric redshifts and to 
perform SED fitting. Grey concentric circles are observed magnitudes in three HST /WFC3 
UVIS, four HST /ACS, three HST /WFC3 IR, one VLT Ks, and two Spitzer /IRAC [3.6], 
[4.5] bands. [Bottom panel] Same as the top panel but for the comparison sample of LBGs 
at z ~ 4-5. Stellar masses are in solar units and stellar ages are in years. 
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Fig. 2. — Comparison of the SED-based photometric redshifts with the publicly available 
spectroscopic redshifts for LBGs at z ~ 1-3. The concentric circles indicate unreliable 
spectroscopic redshifts based on their quality flags. The dotted line shows the 1-to-l relation 
between two redshifts. The histogram shows distribution of (z sp — z p h/l+z sp ), and the la 
uncertainties in this distribution is ~5%. The dot-dash lines show 3cr limits. 
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Fig. 3. — Example fitting of the UV spectral slope (5 to best-fit SEDs for two LBGs from the 
z ~ 1-3 sample. The dotted black line is the best-fit SED, while solid black line shows the 
estimated UV continuum slope (3. The black filled circles are observed magnitudes. We test 
two different rest-frame UV wavelength range (1300-1900 A and 1300-3400 A) to address 
the robustness of our estimated UV continuum slope f3, and find similar median slopes within 
la scatter. 
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Fig. 4. — UV spectral slope j3 as a function of redshift. Blue filled squares indicate the 
median values when we fit rest-frame 1300 to 1900 A wavelengths, and blue open squares 
indicate 0s when we fi t 130 to 3400 A. The higher re dshift (z > 2.5) measurements from 
Finkelstein et ajj (|2012l) andlBouwens et al.l ( 120091 . |2012[ ) are also shown for comparison. Ours 
and the Finkelstein et al.l (120121 ) uncertainties are the standard error of the mean, while un- 
certainties from Bouwens et al. are la scatter. For comparison, our la uncertainties/scatter 
are listed in Table [TJ 
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Fig. 5. — UV spectral slope (3 versus best-fit SED parameters (E(B-V), stellar mass, stellar 
age, star- formation rate). The size of the circle increases with redshift i.e., the smallest circle 
correspond to the lowest redshift bin (z ~ 1.6), and the largest circle correspond to z ~ 4-5 
bin. The changes in the UV-continuum slope (3 from z ~ 1.6 to 5 are most likely due to 
the change (factor of ~2 or 0.3 dex) in the dust content E(B-V) of the galaxies, as other 
parameters show smaller variation as a function of redshift. 
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Fig. 6. — Distribution of stellar masses, stellar ages, star-formation rates, and dust extinction 
for LBGs at z ~ 1-3 (black) and the comparison sample at z ~ 4-5 (red). All four parameters 
are estimated from the best-fit SEDs. Median values are shown by vertical dashed lines and 
la uncertainties are shown by an error bar at the top of the black histogram. Though the 
higher redshift LBGs have, on average, lower values of most of these physical parameters 
than LBGs at z ~ 1-3, the median values are similar within estimated uncertainties. Stellar 
ages do not show any clear evolution with redshift mainly because of higher uncertainties in 
their measurements (see § 13. 3j) . The histograms show normalized numbers for both samples 
(total number of galaxies in each sample is shown in Table H]), and a two-sided K-S test - 
in each panel — indicates a probability P < 0.006 that the two distributions are drawn from 
the same parent distribution. 
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Fig. 7. — Stellar mass versus UV absolute magnitude relation for U225-, U275-, t/336-dropout 
samples. Bottom right panel shows the same relation for the comparison sample of LBGs at 
z ~ 4-5 (black circles for z ~ 4, and black squares for z ~ 5). The average uncertainties are 
shown in the bottom-left corner. The dotted black line is the best-fit line with a logarithmic 
slope of 0.46. The dot-dash line shows the la scatter (~0.3 dex) from the best-fit linear 
relation. 
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Fig. 8. — Stellar mass versus SFR relation for U225-, U275-, L7 336 -dropout samples. Bottom 
right panel shows the same relation for the comparison sample of LBGs at z ~ 4-5 (black 
circles for z ~ 4, and black squares for z ~ 5). The average uncertainties are shown in the 
bottom-right corner. The dotted black line is the best-fit line with a logarithmic slope of 
0.90. The dot-dash line shows the la scatter (~0.6 dex, ~0.4 dex for z ~ 4-5) from the 
best-fit linear relation. In the upper panel, few galaxies form a sharp edge towards high SFR 
values, which could be an artifact due to lower limits on the model parameters r and t (e.g., 
Hainline et aPboid ). 
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Fig. 9. — Distribution of mass-to-light ratios and specific SFRs for LBGs at z ~ 1-3 (black) 
and for the comparison sample at z ~ 4-5 (red) . Median values are shown by vertical dashed 
lines and la uncertainties are shown by an error bar at the top of the black histogram. 
Though the higher redshift LBGs have, on average, lower values of these parameters than 
LBGs at z ~ 1-3, the median values are similar within estimated uncertainties. The his- 
tograms show normalized numbers for both samples (total number of galaxies in each sample 
is shown in Table [TJ, and a two-sided K-S test — in each panel — indicates a probability 
P < 0.05 that the two distributions are drawn from the same parent distribution. 



-29 - 



Table I. UV Spectral Slopes (/?) 



Rcdshift 


Number of 


UV Slope" 


la 


SEM 6 


< z > 


Galaxies (N) 


P 


scatter 




1.6 (±0.2) 


47 


-1.59 (-1.67) 


0.29 (0.27) 


0.04 (0.04) 


2.2 (±0.3) 


126 


-1.71 (-1.81) 


0.34 (0.37) 


0.03 (0.03) 


2.6 (±0.4) 


213 


-1.71 (-1.78) 


0.47 (0.50) 


0.03 (0.03) 


3.8 C (±0.5) 


182 


-1.88 (-1.93) 


0.45 (0.41) 


0.03 (0.03) 



a The UV continuum slope j3 is estimated from best-fit SEDs by fitting 
a line between the rest-frame wavelengths 1300-1900 A (1300-3400 A) 

b The standard error of the mean = (la scatter) / y/N 

c The comparison sample of -B435- and V606-dropouts. Because of small 
number of V(306-dropouts, the average/median redshift is similar to a 
-B435-dropout. 



